A darkless spacetime 
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In cosmology it has become usual to introduce new entities as dark matter and dark energy in 
order to explain otherwise unexplained observational facts. Here, we propose a different approach 
treating spacetime as a continuum endowed with properties similar to the ones of ordinary material 
continua, such as internal viscosity and strain distributions originated by defects in the texture. A 
Lagrangian modeled on the one valid for simple dissipative phenomena in fluids is built and used for 
empty spacetime. The internal "viscosity" is shown to correspond to a four-vector field. The vector 
field is shown to be connected with the displacement vector field induced by a point defect in a four- 
dimensional continuum. Using the known symmetry of the universe, assuming the vector field to be 
divergenceless and solving the corresponding Euler-Lagrange equation, we directly obtain inflation 
and a phase of accelerated expansion of spacetime. The only parameter in the theory is the "strength" 
of the defect. We show that it is possible to fix it in such a way to also quantitatively reproduce the 
acceleration of the universe. We have finally verified that the addition of ordinary matter does not 
change the general behaviour of the model. 

PACS numbers; 98.80.-k, 95.36.+X, 95.30.Sf 



I. INTRODUCTION 



When studying the universe as a whole we have to 
take into account a number of observed behaviours and 
of physical constraints. It is well known that, since its 
early years. General Relativity (GR) provided cosmolog- 
ical solutions able to describe the large scale evolution 
of the universe from a singular event (the Big Bang) to 
the present epoch. While accumulating evidence, how- 
ever, more and more details are entering the scene and 
the original theoretical framework is no more enough to 
account for all of them. This is the reason why people 
have tried and are trying to introduce new theories of 
gravity, alternative to the original GR, or to modify it in 
a way or another. 

Since the moment when growing evidence supported 
the idea that the universe is undergoing an accelerated 
expansion [1] Einstein's "big blunder" [2], the cosmo- 
logical constant, has seen a big revival and has been en- 
riched with new and sophisticated theories. One is led 
to think that the universe is filled up of something ex- 
erting a negative pressure, responsible for the accelera- 
tion. This "something" has been called in various ways, 
the most popular being dark energy which becomes for 
instance quintessence |3] or phantom energy |4] accord- 
ing to some specific theories. The case of the accelerated 
expansion is not the only one to be treated by means 
of some new field. The homogeneity of the cosmic mi- 
crowave background (CMB) radiation seems to imply, 
very close to the big bang, a phase of extremely fast ex- 
pansion and this is accounted for by means of an ad hoc 
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scalar field, the inflaton, with an even more ad hoc po- 
tential [5]. We also know about dark matter, needed to 
explain inhomogeneities of the CMB, the rotation curve 
of galaxies and the behaviour of galaxy clusters [6]. Yet 
another approach consists in using a modified Hilbert- 
Einstein action integral, expressed in terms of some non- 
linear function of the scalar curvature [7, 8]. 

The situation, even though being richer and more var- 
ied, resembles the one with ether at the end of the XIX 
century, and Occam's razor seems to be left aside for a 
while. 

Here we try a different approach taking advantage of 
analogies with other branches of physics. The power 
of analogical deduction has played an important role in 
the past and still proves to be fruitful even today, for in- 
stance in the case of black holes and Hawking radiation 
[9]. 

Our current vision of the cosmos, especially in GR, 
is essentially dualistic, the actors being spacetime on 
one side (left hand side of the Einstein's equations) and 
matter-energy on the other (right hand side of the equa- 
tions). Structures, differences, variety of features be- 
long to matter-energy. The only intrinsic property of 
spacetime, besides the ones induced by matter-energy 
through the coupling constant G, is expressed by the sig- 
nature of the metric tensor. 

The paradigm we are proposing here considers a 
spacetime endowed with some more features on its own 
that remind those of a physical continuum. Whenever, 
in a given physical system, we find a symmetry, we 
know that something real must be there to induce that 
symmetry. In the case of the whole imiverse, its global 
symmetry, in four dimensions, implies the presence of 
a singular event: the center of symmetry. We may state 
it either way: telling that the symmetry implies a zero- 
dimensional singularity, or that the singularity induces 
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the symmetry. The novelty in our approach is in think- 
ing that the singularity is not due to the content (mass- 
energy) of the spacetime, but is built in the very space- 
time. The next step consists in interpreting and treating 
the singularity just as a defect in a continuous medium 
is in the classical elasticity theory Jl^llJ]. A point defect 
in an otherwise homogeneous and isotropic medium in- 
duces a strained state, characterizable by means of a ra- 
dial vector field: the "rate of stretching" in the radial 
position. That vector field is divergence-free, except in 
the center of symmetry. 

In our theory it is indeed a four-vector field that plays 
an important role and it is not at all the first time that 
a vector field is introduced in the description of the be- 
haviour of the cosmos. An example is Bekenstein mod- 
ified theory of gravity (a proposed relativistic version 
of Milgrom's MOND [12]) which can be thought of as 
a scalar-vector-tensor theory [13]. Besides this, the core 
of the so called Einstein-aether theory (i^i-theory) is a 
timelike unit vector field fl^ which is introduced in the 
spacetime action integral in the way this is usually done 
for additional field components. In the Einstein-aether 
theory, in particular, the cosmic vector field appears in 
the Lagrangian through its first covariant derivatives 
and a Lagrangian multiplier related to the unit norm 
constraint. The vector field brings about a violation of 
the Lorentz symmetry, a priori introducing a preferred 
local rest frame, ^-theory has been explored in its con- 
sequences (see for instance 11511 , III61I , and the references 
quoted in [14]) in order to fix upper bounds on its pa- 
rameters and concluding that it would be compatible 
with the present state of observations and experiments 
flSJ. Our approach is however different; for us, the pri- 
mordial fact is a pointlike defect in spacetime and it is 
this defect that induces both the symmetry of the phys- 
ical manifold we call spacetime, and the vector field. In 
this way the form of the field is determined by the very 
existence of the defect, just in the same way as a strain 
state is induced by defects in material continua. 

The idea of spacetime as an elastic continuum with 
properties depending on its "microscopic structure" has 
a story on its own, with illustrious antecedents, such as 
Sakharov IITtII , who tried to explain the "rigidity" and 
Lorentz invariance of the "medium" in terms of zero 
point fluctuations of quantum fields. Presently, Loop 
Quantum Gravity theories [18] consider a sort of atom- 
istic description of spacetime. On the other side, the 
theory of defects is a well established one and was de- 
veloped more than one century ago together with the 
elasticity theory; the idea of extendin g it to more than 
three dimensions is not new lfTll[T9l [20il2lil22ll. but it has 
never been considered more than a curiosity. Recently, 
the idea of some "solid behaviour" has been called in for 
dark energy [23]; there, however, the solid is the dark 
energy itself and its description is the one of a three- 
dimensional medium evolving in time. Our continuum 
instead is the very spacetime and its solid properties are 
described in four dimensions. As it is the case for clas- 



sical GR, we study a global equilibrium state (including 
strain, distortion and whatever else) from the center of 
symmetry (the defect) to infinity. 

In order to write down the appropriate action inte- 
gral for the spacetime considered above (including the 
defect) we remark that the phase space of the system 
is bidimensional, the generalized coordinates being the 
scale factor and its rate of change. A similar phase 
space, whose generalized coordinates are position and 
velocity, is the one describing the motion of a massive 
point particle across a viscous fluid. From this starting 
point, we are able to write down and then generalize 
an appropriate Lagrangian. What we obtain in the end 
is a spacetime displaying inflationary expansion in the 
neighbourhood of the center of symmetry (i.e. the Big 
Bang), then a deceleration-acceleration-deceleration se- 
quence. 

The theory does not exclude the actual presence of 
matter-energy; in order to study its effect on the be- 
haviour of the imiverse as a whole, we introduce it in 
the form of an ordinary perfect fluid, as usual. We find 
that in the negligible pressure era (today) the presence 
of matter has no influence on the global solution. In the 
radiation dominated era the general behaviour is pre- 
served when the matter-energy content is lower than a 
critical value. The theory contains one free parameter, 
which is the "size" or "strength" or "charge" of the sin- 
gularity. We may fix it in such a way that the present 
value of the Hubble constant as well as the age of the 
universe are reproduced. Doing so, we see that the cur- 
rently estimated content of matter-energy in the cosmos 
is well below the critical value, and the position and du- 
ration of the accelerated expansion are consistent with 
the data from observations. 

The paper is organized as follows. In Sec.|ll]we study, 
from the viewpoint of variational methods, the simple 
classical problem of a particle moving in a dissipative 
medium and extend it to the relativistic formalism. In 
Sec. |III]a "dissipative" Lagrangian for spacetime is in- 
troduced, and then specialized to the case of a homo- 
geneous and isotropic empty universe; in Sec. IIVI we 
analyze the effect of the inclusion of ordinary matter; 
Sec. |V] verifies the existence of the Newtonian limit of 
the theory. Finally Sec. I VII is devoted to the summary 
of our findings and to the discussion of our conclusions. 
The signature used in the paper for the metric tensor is 
{+,-,-,-)■ 



II. A CLASSICAL MODEL: A PARTICLE IN A VISCOUS 
MEDIUM 

Once one has decided to try and account for a given 
observation (in our case it is the accelerated expansion of 
the universe) by modifying a consolidated theory, such 
as GR, one needs some criterion to decide what change 
to introduce and test. A possible and often used ap- 
proach is to explore mathematical variants of the ba- 
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da dx 
dz' dt 



a, X 



Figure 1: Phase diagram of a Robertson Walker universe. The 
dash and dot Une corresponds to an inertial expansion (con- 
stant expansion rate); the dotted line is a decelerated expan- 
sion; the dashed line is an accelerated expansion. Observa- 
tions tell us that the actual behaviour of the universe cor- 
responds to both deceleration and acceleration, in different 
epochs. Relabelling the axes with an x instead of a, and adx/dt 
instead of da/dr, the diagramme describes the state of a point 
particle interacting with a surrounding medium. Dash and 
dot is inertial motion (constant speed); dots is motion under 
the action of friction; dash is an accelerated motion driven by 
the medium. 



sic theory, introducing free functions and free parame- 
ters, then adjusting the parameters and the functions so 
to reproduce the observed results and satisfy the phys- 
ical constraints the problem has. This procedure can 
be more or less fortunate, but often poses problems of 
physical interpretation with the newly introduced func- 
tions and parameters. A different way may be, as al- 
ready mentioned in the introduction, to look for corre- 
spondences with other parts of physics. In fact we know 
that in many cases the same set of equations can gov- 
ern apparently unrelated phenomena. This is the case 
for instance of classical field theory and fluidod5mam- 
ics, elastic waves in solids and fluids and electromag- 
netic waves, conservation laws in general, etc. For this 
very reason we decided to look around for some ordi- 
nary situation displaying the same formal properties as 
a Robertson Walker (RW) universe. Indeed, if we accept 
that the universe has the typical symmetry expressed by 
the RW line element, we see that the state of the uni- 
verse is described by the only scale factor a with an evo- 
lution parameter t (cosmic time). The situation may 
be schematized as in fig. |[T]|, where various evolution 
trends are drawn. A simple transliteration (from a to x) 
converts the evolution of the imiverse into the interac- 
tion of a point particle with an isotropic fluid. 

Exactly the same type of phase space is obtained when 



describing the motion of a point particle interacting with 
a surrounding fluid. Now the interesting feature of this 
analogy is that we know how to write the equations 
governing the motion of a point particle of mass m in 
a viscous fluid. Although dissipative, the problem can 
be treated starting from the Lagrangian 



(1) 



A simpler form of this expression was initially intro- 
duced by Caldirola [24], then by others [25] for different 
purposes. The Euler-Lagrange equation deduced from 
Eq. lHJ is 



V .7 

X + 7X + 



0, 



(2) 



where 7 > may be interpreted as the laminar viscos- 
ity coefficient and ?/ is the turbulent viscosity coefficient. 
Actually, Eq. ||2j represents a viscous motion if for tj > 
the motion is progressive; for regressive motion t] has to 
be assumed < 0. In this approach, the properties of the 
fluid and of the interaction are all contained in 7 and rj, 
that are assumed to be constants, which means that the 
fluid is not affected by the motion of the particle through 
it. Eq. ||2ll is however inadequate since it is not invariant 
for reversal of the x-axis. If vq is the initial velocity of 
the particle, the solution of l|2) is 
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Vq ' I 



(3) 



For 



-00 < Co < 



the solution diverges at some 

finite positive time. We now recast the problem in a rel- 
ativistic way. We shall consider a flat spacetime and in- 
troduce the action integral 



-m / e 

I A 



(4) 



where dr' = cdt' , v' = x' and ds = \/l — c'^/c^dr'. The 
exponent in Eq.||4j is assumed to be a true scalar, i.e. the 
scalar product of two four-vectors 

7 = (a,i6,/5,/5) , 
/ = (ri',x',0,0) . 

The reference frame has been fixed so that the x' axis 
coincides with the direction of motion. Using Cartesian 
coordinates, the form of 7 expresses the expected space 
isotropy of the medium. The invariant associated with 
7 is 



and the four-vector has been assumed to be timelike. 
The Lorentz-invariant form of Eq. @ is now 



S = -m r e'?«'''>'"^''ds , 
J A 



(5) 
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The Euler-Lagrange equation from @ is 



(6) 



Everything becomes more transparent and simpler 
if, applying an appropriate Lorentz transformation, we 
change the reference frame so that 



7= (X,0,0,0). 



(7) 



We see that in this case a privileged reference frame ex- 
ists: it is the one of the fluid (unprimed quantities). The 
equation of motion is now 



1 



0. 



(8) 



Equation||8) represents the relativistic version of motion 
in presence of laminar viscosity; now the solution is a 
decelerated motion and no troubles arise from any re- 
versal of the space axes. It is explicitly: 



x = ±- 



\ 



(9) 



+ 1 



M 
e c 



Starting with an initial value Vq < c, the velocity be- 
comes zero in an infinite time; "photons" do not interact 
with the medium (their velocity stays equal to c). We 
could reasonably introduce a dependence of x on v, but 
to discuss further this elementary situation is out of the 
scope of the present paper. What matters is that it is pos- 
sible to give a Lagrangian treatment of a simple dissipa- 
tive phenomenon describing a non-uniform evolution in 
time, in a relativistic context. 



III. THE BEHAVIOUR OF SPACETIME 

We may now use the model in the previous section as 
a guiding idea (by analogy) and the action ||5ll as a sug- 
gestion or inspiration to describe a spacetime that un- 
dergoes expansion or contraction at a non-uniform rate. 
We stress that not the whole universe but the only space- 
time is considered. 

The starting point is the usual Einstein-Hilbert action 



n2 



Rdn, 



(10) 



where R is the scalar curvature and dCl = y^j^d^x is the 
invariant volume element. 

We directly introduce in the spacetime the kind of 
symmetry we usually attribute to the universe, i.e. four- 
dimensional isotropy around a given event. As it is well 



known the most general symmetric line element of this 
type is the RW one 



dT^ - a {tY 



_df_ 



+ r' 



(de^ + sin^ 9d(p' 



(11) 

where k = 0,±1 and r is a dimensionless coordinate 
(product of usual length times the square root of the 
space curvature). Introducing the metric tensor implicit 
in Eq. ((Till i^ito Eq. l ITOl l one has 



-6V, 



ad + a +k] adr , 



(12) 



where Vj: = Jq Jq Jq" -^^^drdOdcp and dots represent 
derivatives with respect to x. 

A. "Dissipative" spacetime 

To take advantage of the simple example described 
in Sec. Ull we may want to reproduce the same logical 
structure while building the action integral. Of course, 
there are important differences to take into account. In 
Sec. In] we had two "actors" entering the scene, the par- 
ticle and the dissipative medium, and an interaction be- 
tween them. The interaction was mediated by a four- 
vector pertaining to the medium and another one de- 
scribing the state (of motion) of the particle. Now the 
"actor" is unique, spacetime itself, and nothing is mov- 
ing across it. However, we may think to the motion of 
the representative point of the state of a hypersurface 
labelled by a in a bidimensional phase space where the 
independent variable is the parameter a. Again, we as- 
sociate to the system a four-vector 7 which couples to 
the state variable in the same way as in the simple ex- 
pression ||5ll. There the position coordinate of the particle 
appeared in the Lagrangian through its (first) derivative 
dx/dt; now, the Lagrangian is written in terms of the 
(second) derivatives of the elements of the metric ten- 
sor The simplest way to couple a vector to the metric 
in order to generate a scalar is by taking the norm of the 
vector itself, so we are led to conjecture the following 
action integral 

S= /"%±^/'"T''')'"Rdn , (13) 

where 7 is meant to represent an internal property of 
spacetime. The correspondence between | [T3l l and ||5| 
may be better seen considering that in ||5| the position 
vector of the particle is made out of the coordinates of 
the particle, whose choice is free, and the final behaviour 
of the system is independent from that choice. In the 
case of | [T3t , i.e. for spacetime, the "coordinates" used to 
describe the state of the system are represented by the 
elements of the metric tensor and again there is a gauge 
freedom in their choice. As written above, the simplest 



5 



scalar we can build from a vector and the metric tensor 
is the one in the exponent of | [T3t . 

Using the metric in the form (TTt , in the (privileged) 
cosmic reference frame, and the four-dimensional ro- 
tation symmetry, 7 necessarily appears in the "radial" 
form of Eq. 0, so that the action | [T5t reads 



e±^' (flfl + fl2 + fc) adr . (14) 



S = -6V;. 



The X = constant case is irrelevant because it corre- 
sponds to simple Minkowski spacetime (the exponential 
may be absorbed in a rescaling of the time coordinate). 
Non-trivial solutions exist only if 



(15) 



and spacetime is not globally homogeneous and 
isotropic. The observation of the universe apparently 
tells us that space is flat, so we further assume k = 0. 
The final effective Lagrangian is 



aa + a \ a . 



The Euler-Lagrange equation for a (t) is 





- ' 1 


(dC\ 


(as) 




\ da 1 



(16) 



, (17) 



or explicitly 

fl (1 ± Ixx'a) 

n1 



(18) 



[2xh'^ ± ± XX") ± 5axx' + 2 = « ■ 



where x' — dx/da. A trivial and non rnterestrng so- 
lution is obtained for a = constant (Minkowski space- 
time). Other solutions however exist. The above equa- 
tion may be reorganized in the form 



-f{a)a , 



(19) 



where 



= Ml±2axx') ■ ^'"^ 

A first integration leads to 

fl = Ae-.r/KR, (21) 

and finally to 



A Jo 



(22) 



B. Choosing ;t; (a) 

The function x{i^) in our case is something which is 
"given" exactly as the global symmetry is. We shall dis- 
cuss this issue in more detail in the next section. What 
matters here is that ;t(fl) is not deducible from a varia- 
tional principle. Should we try and write its field equa- 
tion from the action of Eq. ((14) , the result would trivially 
be = 0, i.e. empty and flat spacetime. A singularity, as 
well as a symmetry, is here an initial condition and not 
a consequence of something else. This means that we 
need criteria and guesses to think about credible forms 
for x{i^)r consistent with the hypotheses. One reason- 
ably simple criterion for a vector field that is expected 
not to spoil the symmetry we assumed for spacetime, is 
to constrain it to be divergence free everywhere except 
in the center of symmetry/ origin of the cosmic times. In 
fact, any event where the divergence of the vector differs 
from zero may be thought of as a "source" and, since we 
need to preserve homogeneity and isotropy, we would 
have to assume a continuous and uniformly distributed 
source. Being, in this paradigm, 7 a concrete quantity 
strictly related to the presence of a singularity, it seems 
more reasonable to have just one pointlike source in the 
origin. There is however also a different, though some- 
how equivalent, reason for choosing 7 to be divergence- 
less, as we shall see further on (Sec lIIICj ), when consid- 
ering the analogy with a defected solid. 

The null-divergence condition is formally written as 

= 7'';,= (/^7'')^^ , (23) 

where the semicolon represents a covariant derivative. 
The solution of this equation is 



(24) 



In practice, the vector field looks like the "electric" field 
of a point charge Q^, but in four dimensions. Inserting 
this result in Eq. | [T9l l, and expressing a in units of Q, we 
obtain 



a _ 36 ± 24fl'^ + fli2 
a ~ 2fl7(fl6;p5) 

As a consequence, we have 



(25) 



(26) 



Choosing the upper signs in | |26t , as well as in ((T4t , the 
expansion rate acquires real values starting from a fi- 
nite non-zero value of a, then it displays a monotonic 
trend which does not correspond to observational data. 
Choosing the lower signs instead, we see that the expan- 
sion rate a has two extrema at 



flM± = V 3 ± V3 



(27) 
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Figure 2: Plot of da/dr as a function of a, in empty spacetime. 
Two extrema delimiting a phase of accelerated expansion are 
clearly seen. 



Figure 3: Scale factor a as a function of cosmic time t; arbitrary 
units. The initial inflationary epoch as well as an accelerated 
expansion between t ~ 0.41 and t ~ 1.3 are clearly visible. 



The corresponding explicit numeric values for the scale 
factors are 



1.08, 



flM+ = 1-68. 



(28) 
(29) 



Fig. |2] shows the behaviour of da/dr as a function of 
the scale factor a. In any case, the asymptotic behaviour 
when fl ^ 00 is fl ^ 0: this is a never-ending expansion. 
From now on, we limit our consideration to the latter 
choice of signs, so, integrating Eq. l|26)l one has 



,1/2 



^5/2 



(30) 



The corresponding behaviour of the scale parameter a 
as a function of the cosmic time t is shown in Fig. |3] 
Close to the origin the negative exponential factor in the 
integral in | |30l l brings about an inflationary phase, that 
cannot be resolved in Fig. |3j followed by a deceleration- 
acceleration-deceleration sequence driving an unlimited 
expansion. In the neighborhood of the origin the scale 
factor a does indeed grow faster than any power of t. 

We can now fix the scale of the expansion. We know 
that in general 



^0 

— =l+z, 
a 



(31) 



where Uq is the present value of the scale factor and z is 
the redshift of light emitted when the scale factor was 
a. From the numerical values ||28ll and ||29t we see that 



the ratio between the scale factor at the end and at the 
beginning of the acceleration epoch is 



Then, using (|3Tt , we can fix 

1+z.- 



2/ 



1.55. 



1.55, 



(32) 



(33) 



where now z, corresponds to the redshift at am- and z^ 
to that at flM+- Looking at the data from the observation 
of high redshift la supernovae |26], we see that (|33t is 
indeed consistent with an initial value of z, ~ 1.6 (or a 
little more) and a final one Zf 0.6. This result is ob- 
tained independently from the value of the integration 
constant Q and in the absence of matter. 

An estimate of Q can be obtained again from | |3H . If 
we let the am- value correspond to Z; = 1.6, we get im- 
mediately, from (|3Tt and (|28t , aq = 2.81, then numeri- 
cally from | |30l l tq = 3.06. Now, if T is the age of the 
universe, it is 



To 



cT 



(34) 



If we constrain T to be not less than 12 billion years (age 
of the oldest globular cluster stars), we conclude that 



Q > 4 X 10^^ m . 



(35) 



This result should also be confirmed considering the 
present value of the Hubble constant Hq. From 
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Figure 4: Point defect originated removing a solid hyper- 
sphere, then squeezing the hole to a point. 



(lower signs) we see that 



Hn 



.3/2 



T = To 



,1/2^ 



pl/{2«^) 



0.21 . 



(36) 



The above result is adimensional. Introducing the ap- 
propriate dimensions it must be 



Hn 



«4 



1.6 X 10~i^s-^ 



or, using the commonly used units. 

Ho ~ 49 km/s x Mpc, 



(37) 



(38) 



Considering the roughness of the model, this is a reason- 
able number for the Hubble constant, which is currently 
estimated to be 65 km/s X Mpc |27]. 

We would like to stress that all this comes from an in- 
ternal property of spacetime, with no matter inside. Mat- 
ter must be further added to the Lagrangian in the tra- 
ditional (additive) way and with a minimal coupling to 
spacetime via the metric tensor. 



C. What could the vector field represent? 

As mentioned above, the "internal" vector field associ- 
ated with empty spacetime can be interpreted by means 
of another analogy with ordinary physics. We know that 
the intrinsic metric of a material continuum can be non- 
Euclidean (non-zero intrinsic curvature) when defects 
are present (see for example lIToll or i28i1 and references 
therein). The corresponding theory has been developed 
many years ago, starting with the formal definition of a 
defect given by V. Volterra |29]. The attempt to extend 
the theory from material elastic media to spacetime has 



been made by many a scientist llTTl[T9ll20ll2lll2"2L I3C 
in various epochs, without leading to a complete formal 
new theory. The similarities are indeed tempting. What 
is easily seen is that, whenever a portion of a continuum 
is removed (or more is added) each point in the mate- 
rial is displaced to a new position (in the unperturbed 
original reference frame) |31i1 



(39) 



The new coordinates are obtained by means of a vector 
displacement field In the continuum a new metric is 
now induced, which is not the original Euclidean one 
but 



g'^=<5'/ + 2£^ 



where 



2 I 3x7 ax' ax/ 



(40) 



(41) 



is the (non-linear) strain tensor. It is important to re- 
mark that the new metric, as well as all physical quan- 
tities of this description, can equally well be expressed 
in terms of the original, undeformed "Lagrangian" coor- 
dinates x„ or of the new intrinsic coordinates x^', being 
the old and the new coordinates numerically identified 
llO'l. In both cases any point is labelled by the same set 
of numbers (the coordinates) plus a vector (the displace- 
ment vector at that point) which is actually zero in the 
unstrained manifold and non-zero in the strained one. 

This framework can be generalized to four dimen- 
sions and to spacetime. The Euclidean basic metric ^ 
5^ is then replaced by the one of Minkowski ri^y and the 
induced metric is written as [31]: 

g}iv = Vfiv + 2£,,v . (42) 

Without further details, let us consider an unperturbed 
(i.e., Euclidean) 4-dimensional space. Then, let us sup- 
pose we remove a 4-sphere and close the hollow by 
pulling radially on each point of the hypersurface of the 
hole. The situation is described in Fig. |4] This proce- 
dure induces a radial displacement field represented by 
a radial four-vector ^. Remarkably, solving the equa- 
tions of the elasticity theory with these symmetry con- 
ditions gives, for ^ in four dimensions, precisely a result 
like Eq. l|24l : the four-vector ^ has a null divergence (see 
for instance Ref. lIToll , Vol. 3 page 107). For spacetime, 
which implies a Wick rotation in order to produce the 
right signature, the induced interval we obtain in these 



^ Here we have used the standard notation for an n-dimensional Eu- 
clidean space with latin indices ranging from 1 to n; on the other 
hand, for 4-dimensional spacetime, the usual greek indices are used, 
going from to 3. 
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conditions, once expressed in appropriate coordinates, 
corresponds to a typical Robertson- Walker metric 



that /" be independent from a ^; in practice it must be 



ds^ = dr^ - (t) dxp^ + ip^ (dO^ + sin^ Odcp' 



, (43) 



where (t) is a non-trivial function and space is flat. In 
order to better explain this result, let us start from the 
general form of the line element of a Minkowski space- 
time expressed in four dimensional polar coordinates 



is = da 



dip^ + ip^ (^de^ + sin^ 0d^2^ , (44) 



where a is now the radial coordinate. The point de- 
fect produces, as written above, a purely radial displace- 
ment field, which means that the only non-vanishing el- 
ement of the strain tensor ||4H is The induced metric, 
according to l|42)l , is then 



dxp^ + ip^ (de^ + sin^ 0d(p^ 



(45) 

Redefining the radial (actually time) coordinate so that 

dT= (^ + da , (46) 

the old radial coordinate a is expressed as a function 
of the new time and the line element becomes | |43l l, as 
claimed. 

In order to clarify the meaning of the 7 vector, we ap- 
ply a procedure typical of the elasticity theory. In a de- 
formed medium the strain is of course accompanied by 
a stress. In linear theory, which we now consider for 
simplicity, the stress tensor a^iy depends linearly on the 
strain tensor e^y (Hooke's law). In our case the chosen 
symmetry implies that the radial-radial components of 
both tensors (i. e. aaa and £„«) be proportional to each 
other. The next step is to think of a given solid angle 
centered at the singularity, then isolate a portion of it 
delimited by two transverse (orthogonal to the radius) 
(hyper)surfaces. In an equilibrium state the forces on 
opposite faces of the boundary of the envisaged piece 
of material must be equal in strength. By definition 
the components of the force on a small surface (four- 
dimensional space) are 



= ^a^^ey^fi^dx^ A dxl^ A dx^ . 



(47) 



Calling again in the symmetry, we see that on the 
"bases" of our portion of solid angle l|47t becomes 



/" = JCe'"'a^dedcj)dip , 



(48) 



where /C is a constant and 6, cp, ip are angles. Now we see 
that the equilibrium within a given solid angle implies 



£""0.^ = constant . 



(49) 



Eq. | |49t is in fact a conservation law. Introducing the 
four-vector 



7 = £ ■ n 



(50) 



where n is a unit vector orthogonal to a given surface 
and 7 represents the flux density of strain. We see that 
the flux of 7 across any closed surface is zero, i.e. V ■ 7 = 
0. Let us identify the 7 in | |50t with the old one, so its 
radial component (the only non-zero component, in our 
case) will be, as before and now by virtue of | |49t and 
(B, 



a^ 

The extended elasticity theory helps us also to find a 
meaning to the Q constant. We know that the energy 
needed to close a void is given by the product of the 
pressure times the squeezed volume W = pV. In our 
case, with the help of the pictorial view of the situation 
shown on Fig. HI we see that the equivalent of the "en- 
ergy" is proportional to 

where A is the radius of the initial hollow. Calling in | |49l l 
and | |50t we see that the "energy" is proportional to (^A. 

is a measure of the ratio between the work done to 
create the defect and its radius. 

This is a consistent logical framework. Part of it relies 
on geometrical bases, whose meaning is clear in space- 
time as well as in three dimensions; this is the case of 
the strain tensor and of the definition (|50l for 7. The 
rest, i.e. the stress tensor with the related quantities, is 
intuitively clear in three dimensions, much less in four, 
however it is a tool for arriving to the final interpreta- 
tion, which remains essentially geometrical. 



D. Equivalent matter distribution 

Once the metric tensor is defined, we can compute 
from it the Einstein tensor and, taking the Einstein equa- 
tions literally, interpret it as being proportional to the 
energy-momentum tensor of some matter-energy distri- 
bution responsible for the peculiar metric. Doing this 



^ When considering the example in the text, the forces on opposite 
sides of the piece of material must of course be opposite in direction, 
but we may consider the force exerted by the external (with respect 
to the singularity) medium on a given surface and in this case the 
direction is everywhere the same. 
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of those theories the equation of state | |56t would come 
from a different Lagrangian than ours; however here the 
comparison is done only at the final stage. What we 
would like to stress is that such an equation of state, if 
sought directly, would appear to be rather artificial and 
indeed it is, if thought as pertaining to an actual "fluid" 
of any sort. 



IV. THE EFFECT OF MATTER 

Let us now verify what the effect of matter is in a 
spacetime like the one described before. To this aim, we 
consider the conceptually simplest situation and intro- 
duce a perfect fluid minimally coupled to the geome- 
try, so that the total Lagrangian of the problem is (see 
Ref . ll32ll and references therein) 



ni 



n, 



pdn , (57) 



Figure 5: Effective pressure as a function of the scale parameter 
a. Arbitrary units. The initial negative values correspond to an 
inflationary epoch. 



exercise in our case produces the following effective 
energy-momentum tensor: 



e^'"' , (51) 





= ^ 


_3 fa 


t 




K 


~ K\a 








laa + 






K 










= T 





kQ2 + 6 

6 (5a'' -6) 
fl3 (^6 + 6) 



(53) 



where it is k = S/rG/c^. 

This energy-momentum tensor has the appearance of 
the one of a perfect fluid whose effective matter-energy 
density is 



P = 



1 



kQ2 + 6 
The corresponding effective pressure is 
6 5a^-6 



(54) 



(55) 



and is represented in Fig.|5] The initially negative values 
correspond to inflation. From Eqs. | |54t and | |55t one can 
immediately obtain the Equation of State (EOS) of this 
fluid 



_ 2{5a^ - 6) 



(56) 



The pressure stays negative up to fl^ = 6/5. In the lan- 
guage of dark energy theories this is equivalent to a pe- 
culiar choice of the factor w = p / p. Of course in the case 



where p is the pressure of the matter-energy fluid. Fol- 
lowing the traditional approach, we consider that in the 
present time the fluid is reduced to an almost incoher- 
ent dust, i. e. p ~ 0. In this condition and with the RW 
symmetry the matter energy density scales as p = po/d^ 
(matter conservation) so that its contribution to the La- 
grangian is simply a constant: the expansion law is un- 
affected. When the presence of the fluid is relevant is 
in the early epochs where the matter density is assumed 
to be negligible with respect to the pressure (radiation 
dominated universe). Conservation of entropy together 
with matter brings about a pressure that scales as 



(58) 



where ipisa positive parameter and k has been included 
for convenience. From l|57|l and in the case of the RW 
symmetry, we obtain the Euler-Lagrange equation: 



36 



+ 1 



24 



. (59) 



Condition (|23l has again been imposed on the 4-vector 
7, as before, so that Ip4) holds. 

Looking for solutions that, in the absence of matter, 
reduce to the already known case | |26t we pose 



a=f{a)e^ 



(60) 



Here f{a) is a function of the expansion parameter a, 
and A is a constant to be determined later Differentiat- 
ing Eq. (|60l with respect to cosmic time t gives 



«=/(/'-6/A).2A/„^ 



(61) 



where /' = df /da. Introducing Eqs. | |6T] | and | |60t into 
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Figure 6: Behaviour of the expansion rate of the universe in 
the presence of ordinary matter: the soHd line shows the de- 
pendence of a on the cosmic time in the case of a spacetime 
with matter in subcritical conditions with \p = 0.5. For the 
sake of comparison we also show the empty spacetime case 
{ip = 0, dashed line) and a case of supercritical matter den- 
sity {tp = 1, thick grey line). We note that the three curves 
start with an accelerated expansion phase. For later times this 
is converted into deceleration, but for \p = 0.5 and ip = an 
effective "re-heating" occurs. 



(imaginary expansion rate). At a smaller scale evidently 
some more refined picture is needed. 

Looking for the extrema and differentiating (|65l one 
obtains the condition: 
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lAa"" +36) + IxpW -18- a^^] =0. (67) 



For )p = the solutions of (|67t are l|27l . Studying the 
equation for ip > we see that three real positive roots 
exist as far as < < ipc, where ipc — 0.8; this means 
that a has three extrema. For ip > ipc only one extremum 
exists. Fig. i flOl l compares the behaviours of an empty 
spacetime with those of one with, respectively, a sub- 
and a super-critical matter content. 

In practice, when matter is present in the form of a 
radiation fluid, the model starts working from a typi- 
cal value a — ip of the scale factor. Initially one has 
a phase of inflationary accelerated expansion, then the 
expansion rate starts decreasing, but after a while, if it 
is ip smaller than the critical value, a sort of re-heating 
happens and the universe accelerates again its expan- 
sion; finally the expansion rate decreases once more un- 
til reaching the value at infinity. When tp > ipc the ini- 
tial accelerated expansion is followed by a never ending 
deceleration. 

The parameter ip scales as Q^, thus, using the estimate 
of Eq. < |35l l, we see that the critical value, in international 
units, is 



Eq. | |59t gives 

6 \ r, 2\la^ 



2 a 



+ 



ffe 



^ ^ 36 - 72A _ 24 + 12A \ _ i^^y.s 



7I2 



Choosing A = 1 / 2, this equation becomes 



Iff [a 



30 



= 



The solution of 



IS 



2 _ Aa^ - xpa^ 
^ 6 + fl6 ' 

A is an integration constant. 
Finally l|60l tells us that the expansion rate is: 



(62) 



0. 



(63) 



(64) 



a = a 



2 f u — ^ 
6 + fl6 



1/2 



(65) 



A comparison with | |26l l fixes A = 1 and the overall sign 
of the formula. 

From l|65)l we see that the model fails to describe the 
situation for 



0<a<xp 



(66) 



10^0 m2 



(68) 



Should we conjecture that the minimal ap value, below 
which the classical fluid description fails, is the Planck 
length, it would be 



ipp ~ 10"^° m^ 



Actually one has that (see for instance BSH ) 



(69) 



(70) 



where pr^ is the present radiation energy density in the 
universe and aq its present scale factor. One usually es- 
timates that jD,-(, ^ 10~^^ J/m^; using for aq the order of 
magnitude of Q we obtain 



!/; ~ 10^^ m^ , 



(71) 



or, in the adimensional form used throughout the paper, 
^ 10^^*^, well inside the subcritical region. The cor- 
responding minimal scale of the universe (below which 
the model is not able to describe what happens) would 
be 



flm ^ lO^'^m. 



(72) 



The conclusion of this section is that the presence of or- 
dinary matter apparently does not spoil the results ob- 
tained in Sec.|III]for empty spacetime. 
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V. THE NEWTONIAN LIMIT 

Of course our theory, as any cosmological theory, 
must prove to be able to reproduce the known results 
at the scale of the Solar system and weak gravitational 
field, which means that it should possess a Newtonian 
limit. 

In order to prove this it is convenient to start from 
the action integral | |57t and write the general form of the 
Euler-Lagrange equations of the theory: 



(73) 



1 -^»,7V^^^ [?/^v] 



2 

kT 



The terms on the second and third line are symmetrized 
in ji and v and T^y is the energy-momentum tensor of 
matter. 

Now suppose that 

where Tfiy is the energy momentum tensor of the cosmic 
fluid and Tfiy is the one of a bunch of matter, we assume 
for simplicity to have stationary space isotropy around 
any given point. 

Under these conditions we expect the perturbed line 
element (spacely isotropic coordinates) to be: 



ds^ = (1 + ho) dT^ - fl^ (1 + hs) (dx^ + dy^ + dz^ 



(74) 



with ho, hs << 1 and depending on r = y/ x"^ + y"^ + z"^ 
only. 

The perturbed metric tensor will perturb the flow 
lines of the vector 7 field also, inducing the same kind 
of space symmetry. So we write for the components of 
the perturbed vector: 



Y° = x(l+/o) 
x' 



(75) 



where the f's are assumed to depend on r only (as the 
h's from which they stem) and fo, /s < < 1 (at least as 
small as h's). 

Recalling that the Y field has its origin in the cosmic 
point defect and that no other defect has been intro- 
duced, the divergencelessness condition Y?J, = must 
still hold. Using the equivalent form (-^— gV') = 0, 

the zero order (unperturbed) solution x = 1 /fl^/ l(74|l . 



and the dependences of the f's and h's, we obtain in the 
first order approximation 



which implies 



/« = 3 







(76) 



(77) 



b is an integration constant. 

A further constraint we can introduce is that the norm 
of the Y vector remains unchanged. This is because, 
again, the vector field depends solely on the existence 
of a defect and the global symmetry it induces. Consid- 
ering this constraint we write 

Yf,Y'' = 7,,7?' = 
or, using (|74l , (|75t , and stopping at the first order. 



which implies 



X^ = X^ {l + 2fo + ho) 



fo 



(78) 



Once these constraints have been implemented we 
may go back to l(73l and consider the time-time equa- 
tion: 



e-^' (Goo - jIr) 



- gooe-^' (27>.;<.7A7^'" " 7''^"7«;. " I'l^c. 



''gaO (27"7«;07A7^'" - 7"'"7«;0 - 7"7«;0'" 



kToo. 



The metric tensor is diagonal, which fact implies 
e-x' (Goo - 7oR) 



kT, 



00- 



The next step is to expand everything up to the first or- 
der in h's and f's; drop the zero order terms, which are 
satisfied by l|24t and | |65t with the source Tqq', use condi- 
tions ((78|l and l(76|l . The remaining first order equation 
is: 

*\ 



3%ho 



~l~ 2 i,^s,xx ~l~ hs,yy + ^s,zz) 



+ {h0,XX + ho,yy + /lO,zz) 

~l~ 2^^ {hs,xx ~\~ hs,yy + ^s,zz) 



(79) 
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Now, for ordinary time scales the rate of change of a, 
i.e. the Hubble constant, is extremely small so that we 
neglect the first term in l(79)l . Passing to an orthonormal 
base (marked by a '~) the factors are absorbed into 
the space derivatives, so that l(79l becomes: 

+ (v^/zo + IV^hs) = Ke^Voo (80) 

Finally, exploiting the gauge freedom in the choice of the 
coordinates (Lorentz gauge with time independent /z's), 
(|80l is reduced to 

which can be read as the Poisson equation for a New- 
tonian gravitational potential with a renormalized cou- 
pling constant 



slowly changing in cosmic times. The vector field 7 
appeares at the local scale only through its norm 
included in the renormalization factor of the Newton 
gravitational constant G. 

VI. CONCLUSION AND DISCUSSION 

We have applied a heuristic approach to the problem 
of describing the behaviour of the universe in its expan- 
sion. Instead of introducing new components in what 
should correctly be called "matter" (any scalar or tensor 
field usually considered is indeed "matter" in the sense 
that it contributes to the right hand side of the Einstein 
equations and appears additively in the Lagrangian), we 
have used a model based on the idea that the very space- 
time is endowed with a property analogous to the inter- 
nal viscosity of a fluid. This feature has been treated in- 
troducing an exponential factor in the Langrangian and 
exploiting from the very beginning the four-symmetry 
we think the universe has around the origin. The scalar 
in the exponent of the new factor is thus built from a ra- 
dial (in four dimensions) vector field. Symmetry consid- 
erations suggest that the vector field can be divergence- 
less everywhere except in the origin. Solving the vac- 
uum Einstein equations in these conditions we end up 
with a global RW metric with a scale factor depending 
on cosmic time in such a way to reproduce an initial in- 
flationary era, followed by a decelerated, then again ac- 
celerated, and finally decelerated expansion of the uni- 
verse, or, to say better, of space itself. 

Looking for a possible explanation of the "friction" 
described by the new vector field, we have had re- 
course to a further analogy with ordinary material con- 
tinua. We have assimilated empty spacetime to a four- 
dimensional continuum containing a pointlike defect. 



and then we have analyzed the strain and consequent 
metric tensor induced by such a defect. The final re- 
sult is again a RW metric with a time dependent scale 
factor The radial displacement field of this scenario is 
divergence-free as the strain flux density is, thus allow- 
ing the identification of the latter with the initial "vis- 
cous" vector field. The model described here presents 
results that scale as the"strength" or "charge" of the 
center of symmetry, which may be fixed in order to re- 
produce the expected age of the universe; this param- 
eter is geometrically interpreted as representing both 
the size of the original void and the "rigidity" of space- 
time, where the final defect (initial singularity) comes 
from. Besides this fact, the model, without any further 
recourse to free parameters, reproduces reasonably well 
the observed duration of the acceleration era. Further- 
more, it gives also rise to an initial inflationary era. 

One can wonder what changes this theory brings 
about in the Newtonian limit and on a local scale. It 
is easy to see that with respect to this nothing hap- 
pens. The relevant quantity in the action l(T4l is the 
scalar x^ whose rate of change with cosmic time is 
"^Xf X = — 6fl/fl = — 6H; in practice, with the present 
value of the Hubble constant, one has 2x/ X ~ —10^^"^ 
s^^ whose inverse corresponds approximately to 3 bil- 
lion years. For time intervals much smaller than the time 
scale T ^lO^"^ s the exponential factor in l(T4l is practi- 
cally constant, thus the known results of GR hold. If for 
instance we further introduce in the Lagrangian the typ- 
ical space symmetries of the Schwarzschild problem we 
obtain the corresponding solution with its Newtonian 
limit. Only for time intervals comparable with T one 
can expect changes, which would show up, still using 
Schwarzschild as an example, in the form of an adiabatic 
change in the unique parameter not fixed by the space 
symmetry, i.e. the mass of the source. Many would pre- 
fer to state it as a time dependence of the effective grav- 
itational constant over cosmic times, but the result is the 
same. Summing up, we see that GR appears as a short- 
time approximation of the theory we propose. 

Our final step has been to verify that the addition of 
ordinary matter in the form of a fluid (with the densities 
we obtain from observational data) does not subvert the 
behaviour of the universe we obtained for empty space- 
time. Simply, the compound model (spacetime plus 
matter) starts working from a minimum scale factor, be- 
tween the Planck era and the present epoch. 

The form initially chosen for the action integral is 
somehow reminiscent of other approaches, from string 
theory to f (R) theories, without however coinciding 
with any of them. We ourselves showed how the effects 
may be thought of as being due to an effective fluid with 
a peculiar equation of state. As a matter of fact, we ob- 
tained our result following analogies coming from facts 
of known classical physics and introducing reasonable 
(to us) hypotheses, rather than new ad hoc entities. 

Besides our initial motivation to look for an expla- 
nation of the accelerated cosmic expansion (since our 
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theory is a modification of standard GR), we obvi- 
ously would like to verify what the consequences of the 
new spacetime Lagrangian are not only for the Newto- 
nian limit discussed above, but also for such phenom- 
ena as the propagation of metric perturbations (gravi- 
tational waves), propagation of electromagnetic waves 
(modified Ernstein-Maxwell equations), exact solutions 
in various symmetry conditions (the equivalent of the 
Schwarzschild and Kerr solutions) etc. Since ours is 
a metric theory having Minkowski both as the tangent 
and the asymptotic spacetime we do not expect, at least 
on not too big scales, relevant changes with respect to 
the standard theory. However, the differences could 
show up both in local high curvature regions of space- 



time and on the large scale behaviour of matter systems. 
To explore all these possibilities is our programme for 
the near future. 
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